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Abstract 

An analytical approach, based on the unitary transformation method, has been 
developed to study the effect of quantum lattice fluctuations on the ground state of 
a model electron-phonon system. To study nonadiabatic case, the Green's function 
method is used to implement the perturbation treatment. The phase diagram and the 
density of states of fermions are obtained. We show that when electron-phonon coupling 
constant o? jK decreases or phonon frequency increases the lattice dimerization and 
the gap in the fermion spectrum decrease gradually. At some critical value the sys- 
tem becomes gapless and the lattice dimerization disappears. The inverse-square-root 
singularity of the density of states at the gap edge in the adiabatic case disappears 
because of the nonadiabatic effect, which is consistent with the measurement of optical 
conductivity in quasi-one-dimensional systems. 
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1 Introduction 



The physical and chemical properties of materials with a quasi-one- dimensional charge- 
density-wave (CDW) state, for example, the halogen-bridged mixed-valence transition-metal 
complexes and the conducting polymers, have been the subject of intense study in recent 
years, because of their intrinsically interesting properties and important technological appli- 
cations m, 0. Among the models for one-dimensional systems the Holstein modeip] and Su- 
Schrieffer-Heeger (SSH) model are the two typical electron-phonon coupling Hamiltonian 
studied by many previous authors. The Holstein model is for the on-site electron-phonon 
interaction with dispersionless phonons coupled with electron density operator, while the 
SSH model is for the on-bond electron-phonon interaction. 

In this work we deal with another one-dimensional electron-phonon coupling model, 
the simplified one band Hamiltonian for the halogen-bridged mixed-valence transition-metal 
complexes §|, 0, 



Ls 



where c\ ^ and q^^ are the creation and annihilation operators of electrons at site / with spin 
s, M/(with conjugated momentum pi) is the displacement of the / ion, to is the supertransfer 
energy of electrons between neighboring two orbitals, a is the electron-phonon coupling 
constant, K is the elastic constant and M the mass of ions (throughout this paper we set 
h = ks = 1). We will show that, at least for the spinless case, this model is the same as the 
Holstein model in the adiabatic limit {M oo) but may be similar to the SSH model in 
the antiadiabatic limit (M 0). 

Within the adiabatic mean-field approximation, that is, treating the phonon degrees of 
freedom classically, the model can be solved easily. In the half-filled-band case, the system 
undergoes a Peierls instability and the ground state is dimerized with an energy gap at the 
Fermi points k = ±7r/2||^. However, this approach is questionable and it has been shown 
that the quantum lattice fluctuations must be taken into account to satisfactorily describe 
some physical properties of quasi-one-dimensional systemsP]. Generally speaking, the quan- 



tum lattice fluctuations may decrease the CDW order parameter |[10|, 11]. As the density of 
states (DOS) is concerned, the results of adiabatic approximation have inverse-square-root 
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singularity at the gap edge. By considering the quantum lattice fluctuations, the singularity 
at the gap edge may disappear |jl2|. The relationship between Peierls distortion and phonon 



frequency in the range from cu^ = to cu^ — > oo should be studied for understanding the 
physical properties of electron-phonon interactions in nonadiabatic case. 

In Hamiltonian (|l]) the operators of lattice modes, ui and pi, can be expanded by the 
phonon creation and annihilation operators b^_q and bq, and after Fourier transformation the 
H becomes 

H = + J) + I]efc4,sCfc,s + -yff9{q){bq + &-g)4+g,«Cfc,., (2) 

q ^ k,s q,k,s ViV 

where = —2to cos k is the bare band function, is the total number of sites. The disper- 



sionless phonon frequency = yK/M and the coupling function g{q) = ayl/{2MujT,)[l — 
exp(ig)]. 

2 Effective Hamiltonian 

In order to take into account the fermion-phonon correlation, an unitary transformation is 
applied to H, 

H' = exp{S)Hexp{-S), (3) 

where the generator 5* is 



N q,k,s 

Here we introduce a function S{k',k) which is a function of the energies of the incoming 
and outgoing fermions in the fermion-phonon scattering process. We divide the original 
Hamiltonian into H = + H^, where contains the first two terms of H and the last 
term. Then the unitary transformation can proceed order by order, 

H' = H'' + H' + [S, H''] + [S, H'] + ^[S, [S, //°]] + 0{a'). (5) 



The first-order terms in H' are 

H' + [S,H''] = -^Y.9{q){blq + bq)cl^q^,c,,^ 



q,k,. 
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k+q,s 



^k,s- 



(6) 



Note that the ground state Iqq) of the non-interacting system, is a direct product of a 
filled fermi-sea \FS ) and a phonon vacuum state 



|^o) = |i^5)|AO) 



Applying the first-order terms on \gQ) we get 



{H' + [S,H'])\g,) 



J2 9{q)b-gcl^q,,Ck 



q,k,s 



l-6{k + q,k){l 



(7) 



1^70 ) (8) 



since bq\ph, ) = 0. As the band is half- filled the Fermi energy ep = 0. Thus cl_^_^Ck\FS ) ^ 
only if efc+q > and < 0. So, we have 



{H' + [S,H'])\go) = 0, 



(9) 



if we choose 

5ik + q,k) = 1/(1 + |efc+, - ek\/uj^). (10) 

This is nothing but making the matrix element of + [S, if°]between l^fQ ) and the lowest- 
lying excited states vanishing. Thus the first-order terms, which are not exactly canceled 
after the transformation, are related to the higher-lying excited states and should be irrele- 



vant under renormalization 13 .The second-order terms in H' can be collected as follows: 



[S,H'] + -[S, [S,H^] 



= ^ E E ^-^^^S{k + q, k)[2 - 6{k' + q', k')]{bl^, + - bq) 

q,k,sq',k' 

^ (cI+g^sCfc',s5fc,A:'+g' — cl,_^_g,^gCk,sSk',k+q) 

+ E E + q, mk' + g', k'){e,. - e,,^,,)ibl^, - 6,0(&^ " bq 

q,k,sq',k' t 
^ {ci+q,s'^k',sSk,k'+q' " cl,_^_ , Ck,sSk',k+q) 



1 j2 E ^-^^^^6{k + g, k) [2 - 5(A:' - g, fc')]cU,,.c.,.cI,_,,,,c.^ 



2A^ 



(11) 

q,k,sk',s' 

Sk',k+q is the Kronecker 6 symbol. All terms of higher order than will be omitted in the 
following treatment. 



For the dimerized state, the neighboring atoms move in opposite directions. To take into 



account the static phonon-staggered order ing[|lj], we make a displacement transformation 
to H' 

H = exp{R)H'exp{-R). (12) 

Here 

R=-u-^y^o^l^ibi-k), (13) 

and exp(i?) is a displacement operator: 



exp(i?)«,exp(-/?) = (-l)'uo + Ey^]^^^(^-9 + Mexp(igO- (14) 
Applying the transformation on phonon annihilation operator, we get 



exp{R)bg exp(-i?) = bg + UQ\j ^^^'^ 5^,q- (15) 

If the ground state of if is l^? ), then the ground state oi H \s\g' ): \g) = exp(— 5*) eyip{—R)\g' ). 
We assume that for \g' ) the fermions and phonons can be decoupled: \g' ) ~ |/e ) ), 
where |/e ) is the ground state for fermions. After averaging H over the phonon vacuum 
state we get an effective Hamiltonian for the fermions. 



H,ff = {ph,0\H\ph,0] 



2 

k,s k>0,s 



^ E E ^-^^^^Sik + q, k)[2 - 5{k' - g, fc')]cU.c,,.cI,_,,,,c., (16) 



^ q,k,sk',s' ^' 



where 



Eo{k) = 6,-4E ^^^^ '''h {k\k)6{k,k'){e,-e,,), (17) 



Ao{k) = 2auo[l- 6{k-n,k)]. (18) 

We find by means of the variational principle that the dimerized lattice displacement order- 
ing parameter is 



Uq 



^ E 2«[1 - - vr, k)] Ife (4_,,,Cfc,, + c^c^..,,) fe) . (19) 
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Note that in the adiabatic hmit where a;^ = one has S{k', k) = and H^ff goes back 
to the adiabatic mean-field Hamihonian, 

1 

Kl\Uo + 

k,s k>0,s 



Heffi^Jn = 0) = -KNul + J2 efc4,.Cfc,s + ^ 2auo{cl_^^,Ck,s + c^Cfc_^,,). (20) 



This is the same as the adiabatic mean-field Hamiltonian for the dimerized Holstein model 13 



On the other hand, in the antiadiabatic limit where uj^ oo, we have uq = 0, 6{k', k) = 1 , 
and Heff becomes 

rj f ^ t ^ 9{(l)9{-<l) t t /o1^ 

HeffiuJ^ ^ oo) = 2^ efc4,.Cfc,. - T7 2^ ■ 4+g,.Cfc,.4._g,,/Cfc/,,/. (21) 

k,s k,k' g,s,s' 

Returning to the real space, this Hamiltonian is 

l,s l,s,s' l,s,s' 

(22) 

For the spinless case, the on-site interaction disappears because of the Pauli principle and 
this is the antiferromagnetic XXZ model (through Jordan- Wigner transformation) ||TI| . The 
situation is very similar to the spinless SSH mode in the large u.,^ limit [Q. XXZ model can 
be solved exactly and there exists a transition point at /2K = to[|3- ^'^^ spin-i case, 
(|2^ ) is the negative-t/ extended Hubbard model with V = —U. 

Our effective Hamiltonian works well in the uJt^ = and 00.,^ ^ oo limits and, furthermore, 
for the spinless case He// can be solved exactly in both limits. In this work we concentrate 
on the nonadiabatic effect due to finite phonon frequency uj.,^ in the spinless case because we 
have the exact solutions in the two limits. The last term in H^ff is a four-fermion interaction. 
As we are dealing with a one-dimensional system, how to treat the four-fermion interaction 
is a difficult problem. Since the case for the small u^yr, cUtt < 2to, is very different from that 
for the large c^jr, > 2to; we treat -f^e// in these two cases with different methods. 



3 UJ-jj- <C 2/^0 

In this case the four-fermion term goes to zero as lOtt 0, so it can be treated as a pertur- 
bation and the unperturbed Hamiltonian is 

H'^ff = lKulN + Y.Eo{k)clc, + Y.Aoik){cl_^c, + cW^^). (23) 

^ k k>0 
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The four-fermion term can be re-written as 



1^ ^ 9M9LAsik + q,k)[2-Sik'-q,k')] 
g>0, fc>0,fc'>0 

g>0, fc>0,fc'>0 

q>0, k>0,k'>0 



(24) 



We can distinguish between different physical processes. The first term in Eq.(]2^) is the 
forward scattering one, the second is the backward scattering one, and the last is the Umk- 
lapp scattering one. We use the Green's function method to implement the perturbation 
treatment and it is more convenient to work within a two-component representation, 

/ 




A; > 0. 



(25) 



Thus we have 



(26) 



and the Hamiltonian becomes 



k>0 



k>0 



(27) 



H' 



eff 



-^E E ^-^^^^S{k + q,k)[2~6{k'-q,k')]^l^g^,^l_g^,, 

q>0, k>0,k'>0 
q>0, k>0,k'>0 

1 ^ ^ ^(i^^Ml^,(, + ,,,_.)[2_,(fc'_,,fc'-.)] 



2N 



q>0, k>0,k'>0 



(2^ 
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(^piP = X, y, z) is the Pauli matrix. The matrix Green's function is defined as (the tempera- 
ture Green's function is used and at the end let T — 0) 

G{k,T) = - < T.vl>,(r)vl/t (0) > 

= ^'^(iW{-^^^nr)G{k,uJn)- (29) 

P n 

The Dyson equation is 

G{k,Un) = Go{k,LJn) + Go{k,UJn)T.*{k,UJn)G{k,Un), (30) 

where 

Go{k,uJn) = {iujn - Eo{k)a;, - Ao{k)a^y'^ , (31) 

is the unperturbed Green's function. The self-energy S*(fc,u;„) can be calculated by the 
perturbation theory, 

S*(fc,-n) = ^^^ ^(^^-^)^(^-^\ (A:-,fc)[2-^(A;,fcO] 

X {GQ{k', Urn) + Tr[(7:,Go{k\ t^.m)](Tj 

- mI^I^ 6{k, k - 7r)[2 - S{k, k -7c)] 

X [icryGo{k', uJm)i<yy - a^Goik', uJm)(Tx] 

+ ^ E E ^^^^^^^[S{k, k-7r)+ 6{k', k'-7T)- 6{k, k - 7T)6{k', k' - tt)] 

X {Tr[iayGQ{k',Urn)]i(yy+Tr[OxGQ{k\Urn)]cTx] ■ (32) 

In the perturbation calculation we have taken into account the fact that the forward and 
backward scattering terms contribute nothing to the "charge" gap [|I^, T^. From Eq.( 



one can get that S*(/c,ct;„) is irrelative to c<j„, therefore the spectrum structure of G{k,ujn) 
should be 

G{k,un) = {iu;n-E{k)a,-A{k)axy\ (33) 
From G{k, Un) the fermionic spectrum in the gapped state can be derived 



W{k) = JE^{k)+A^{k). (34) 
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The renormalized band function and the gap function are 

E{k) = i5;o(A:)-|^E W^^V(^^^)[2-^(^^- 



cos^ 



^' Sik' - 71, k)[2 - 5{k' - vr, k)A , (35) 



A(A;) = 2aMo[c-rf5(A;-7r,A;)]. (36) 



Where 



c = 1 + ^E^(^-^,^) , (37) 



™.>o '2aujEl{k) + /\l{k) 



The equation to determine Mq is 

(39) 

In nonadiabatic case Uq is a variational parameter and cannot be measured by experiment 
or Monte Carlo simulation. The quantity which can be measured is the phonon- staggered 
ordering parameter mp, 

= ^E(-i)'<^'> 



"0 - E -^(^ - ^' ^) < 4-7rCfe > 

^E < /e|^Ia,*fc|/e > 
-"■^^ fc>0 



These are basic equations for the u^, < 2to case. If u-,, = we have S{k', k) = and c = 1, 
Eq.(^) becomes the same as that in the adiabatic theory. 

Fig.l shows the phonon-staggered ordering parameter arrip/tQ as function of the electron- 
phonon coupling constant / Kto in the cases of oj-j^/Iq = 0.3 and 0.5. As shown in the 
figure, the dimerization parameter amp/ to increases as the electron-phonon coupling constant 
o? jKto increases but decreases as the phonon frequency ciJ^/to increases. Fig. 2 shows the 



normalized phonon-staggered ordering parameter rrip/mpQ {rripo is the adiabatic value when 
= 0) as functions of the normalized phonon frequency uOt^Iuj-^c iy^-Kc is the value when 
rrip = 0) in the cases of / Kt^ = 0.5 and 1.0. As shown in the figure, the dimerization 
parameter rrip decreases as the phonon frequency cu^ increases but increases as the electron- 
phonon coupling constant a'^/Kto increases, which indicates that the nonadiabatic effect is 
to reduce the dimerization ordering parameters. 

Comparing Eq.(|36|) with that in the adiabatic case, A(A;) = 2auo, we have the gap in the 
nonadiabatic case, 

A = A{n/2) = 2auo[c- d]. (41) 

This is the true gap in the fermionic spectrum. 

From Eq.(l39|), let uq = 0, we get the self-consistent equation of phase transition points 
in the / K ~ uo^^ plane. 

The density of states(DOS) of fermions is 

1 / d 



2vr\rf/tv^'(^) + ^'(^)^=/H) ' (43) 



where, k = f{uo) is the inverse function of = \J E'^{k) + A^(/c). Fig.3 shows the density 
of states(DOS) of fermions for some values of phonon frequency cj^. One can see that a 
nonzero DOS starts from the gap edge and, for small values of uJt^, there is a peak with a 
significant tail below it. The inverse-square-root singularity at the gap edge in the adiabatic 
case|0 disappears. 

4 Ct^TT > 2/;o 

In this case i^e// can be re- written as 

Heff = \KulN + Y.Eo{k)clck + Y.^o{k){cl_^Ck + clck-.) 

^ k k>0 

_ 1 ^ 9_(MzSls(k + q, k)[2 - S{k' - q, fcO]cU,c.cL^c.,. (44) 

q,k,k' 
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According to Eg. (p2D , when u;^ oo, the effective Hamiltonian includes only electron- 
electron interactions of on-site and nearest-neighbor. Therefore, we treat the terms of on-site 
and nearest-neighbor interactions as the unperturbed Hamiltonian and the others as pertur- 
bation because they go to zero when oo. In real space, the unperturbed Hamiltonian 
is 

H'^jf = ^irM^Ar-to[l-i^K)]E(4cm + cJ+iQ) + l^iEcW+iQ+i, (45) 



where 



V, = 2 ^ \9i^)\"^^'^ s^k + q,k)[2-5ik'-q,k% (46) 



^, , ujj, 2o? r'^ , (1 — cos /c cos A;')(cos A; — cos /c') cos 

Fiu^) = / dk / dk- . (47) 

^ (27r)2 K J-n J-n {uj^ + 2to\ COS k - COS k'\)^ ^ ' 



^ ) is the antiferromagnetic XXZ model and can be solved exactly The result of 



Yang- Yang shows that there exists a transition point atP| 



V^i = 2t,[l-F{uj^)]. (48) 

This equation determines the diagram of phase transition in the case of cu^ > 2to- Fig. 4 
shows the phase diagram. We use c<j^/(c<j^ + 2to), instead of cu^, as the variable because it 
goes to 1 when uJt^ oo. The solid line is the result of Eq.(|^) and the dashed line is the 
result obtained from Eq. (^8]) . The dashed-dotted line is the result of equation 



(49) 



K \Lj^ + 2to, 

One can see that, although the formula is very simple, the interpolated result is, at least, 
qualitatively correct. 



5 Conclusions 

The effects of quantum lattice fluctuations on the ground state of a model electron-phonon 
system are studied through an analytical approach. Our results show that when the electron- 
phonon couphng constant / K decreases or the phonon frequency tu^ increases the lattice 
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dimerization and the gap in the fermion spectrum decrease gradually. At some critical value, 
the system becomes gaplcss and the lattice dimerization disappears. A phase diagram in 
the o? I K r-^ lOt^ plane and the density of states of fermions are derived. The inverse-square- 
root singularity of DOS at the gap edge in the adiabatic case disappears because of the 
nonadiabatic effect, which is consistent with the measurement of optical conductivity in 
quasi-one-dimensional systems. 
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Figure Caption 

Fig.l The dimerization parameter arUp/tQ as function of the electron-phonon couphng con- 
stant / = a'^/Kto for cj^/to = 0.3 and 0.5. 

Fig. 2 The normahzed dimerization parameter mp/nipo as functions of the normahzed phonon 
frequency Ut^/ut^c in the cases of / = 0.5 and 1.0. 

Fig.3 The density of states(DOS) of fermions for a'^/Kto = 0.4 with u^/to = 0.001 and 0.01. 
Fig.4 The phase diagram from Eq.(|2D, Eq.® and Eq.®. 
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